plot.extracted <- function(ext.obj, mfrows = c(3,3), wh.vars = NULL, pretty.names = NULL){
  obj.coef <- ext.obj[["coefs"]]
  obj.se <- ext.obj[["ses"]]
  if(is.null(wh.vars)){
    wh.vars <- 2:length(obj.coef)
  }
  par(mfrow = mfrows)
  for(i in 1:length(wh.vars)){
    if(!is.null(pretty.names)){
      main.name <- pretty.names[i]
    }else{
      main.name <- wh.vars[i]
    }
    plot(density(obj.coef[[wh.vars[i]]], na.rm = TRUE), main = paste(main.name, ", SE = (", round(min(obj.se[[wh.vars[i]]]), 2), ", ", round(max(obj.se[[wh.vars[i]]]), 2), ")", sep="")
         )
  }
}
